The divergence of neighboring magnetic field lines and fast-particle diffusion in strong 
magnetohydrodynamic turbulence, with application to thermal conduction in galaxy clusters 
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Using direct numerical simulations, we calculate the rate of divergence of neighboring magnetic field lines 
in different types of strong magnetohydrodynamic turbulence. In the static-magnetic-field approximation, our 
results imply that tangled magnetic fields in galaxy clusters reduce the electron diffusion coefficient and thermal 
conductivity by a factor of ~ 5 — 10 relative to their values in a non-magnetized plasma. 
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The diffusion of fast particles in turbulent magnetized plas- 
mas is important for fusion experiments, cosmic-ray propaga- 
tion, and thermal conduction in galaxy-cluster plasmas (jj, |2j, 
SHEHHUHESmElEi-We consider particle diffu- 
sion in a static magnetic field, which is a reasonable first ap- 
proximation for particles moving much faster than the E x B 
velocity of field lines. The effects of turbulent bulk motions 
and field evolution have recently been considered by 00. 

We consider magnetic fluctuations that possess an inertial 
range extending from an outer scale Ib to a much smaller in- 
ner scale Id with the magnetic energy dominated by scales 
~ Ib- Except where specified, the discussion focuses on the 
case in which there is no mean magnetic field pervading the 
plasma. If a particle is tied to a single field line and travels a 
distance I ^> Ib along the static magnetic field, it takes ~ l/lg 
random-walk steps of length Ib, resulting in a mean-square 
three-dimensional displacement of 



((Ax) 2 ) =al B l, 



(1) 



where a is a constant of order unity. The values of a for the 
turbulence simulations used in this paper are listed in table H] 
When there is a mean field Bo comparable to the rms field, 
((Ax) 2 ) in equation Q is interpreted as the mean-square dis- 
placement perpendicular to Bo. If the particle's motion along 
the field is diffusive with diffusion coefficient D\\, then 



(2) 



and 1151 M 



((Ax) 2 ) ~ t x >\ 



(3) 



indicating subdiffusion [D = lim,^<x,((Ax) 2 )/6? — > 0]. 

In fact, a particle is not tied to a single field line. After trav- 
eling a short distance along the field, field gradients and scat- 
tering cause a particle to take a step comparable to its gyrora- 
dius p across the magnetic field, from its initial field line, Fl, 
to a new field line, F2. If the particle were to follow F2, it 
would separate from Fl because neighboring field lines tend 
to diverge. We call z s the distance the particle would have to 



follow F2 before separating from Fl by a distance Ib- (A par- 
ticle typically separates from its initial field line by a distance 
Ib after traveling a distance slightly less than z s since it drifts 
across the field continuously, but we ignore this effect in this 
paper.) 

If a particle "takes a random step" of length mz s along the 
magnetic field, where m is a constant of order a few, reverses 
direction, and then takes another random step of length mz s 
back along the field, it doesn't return to its initial point. Part of 
the second step retraces part of the first step, but the remainder 
is uncorrected from the first step. This loss of correlation 
leads to a Markovian random walk When mz s ^> fa, a 

single random step corresponds to a 3D displacement of 



(Ax) 2 ~ amzsh- 



(4) 



When mz s ^> A., where X is the Coulomb mean free path, a 
single step takes a time 



A/. 
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When mz s is only moderately greater than Is or X, equa- 
tions @ and @ remain approximately valid. During suc- 
cessive random steps, a particle will find itself in regions of 
differing magnetic shear, and thus z s will vary. The diffusion 
coefficient is given by D = ((Ax) 2 ) /6(At) where (...) is an 
average over a large number of steps 11711 . Ignoring factors of 
order unity, 



D 



"is 



as in 



14, 151, with 



Ls = 7~T- 



(6) 



(7) 



If there is a mean magnetic field comparable to the fluctuating 
field, equation is recovered provided D is replaced by D±, 
the coefficient of diffusion perpendicular to the mean field. 

We treat field-line separation in a static field using 
magnetic-field data obtained from four direct numerical sim- 
ulations of incompressible MHD turbulence. Key simulation 
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parameters are given in table H] The numerical method is de- 
scribed in 1 18]. Each simulation uses Newtonian viscosity v 
and resistivity r| and is run until a statistical steady state is 
reached. In simulations Al and A2, the mean magnetic field is 
zero, and the initial magnetic field is dominated by large-scale 
fluctuations containing 10% of the maximum possible mag- 
netic helicity. Turbulent fluctuations are sustained by non- 
helical random forcing of the velocity field, and the magnetic 
Prandtl number P m = v/r| equals 1. These simulations reach a 
statistical steady state with Kolmogorov-like kinetic and mag- 
netic power spectra. In simulation B, P m = 75 and the field is 
amplified from an initially weak seed magnetic field by turbu- 
lent velocities sustained by non-helical random forcing. When 
the dynamo growth of the magnetic field saturates, the mag- 
netic field is dominated by fluctuations on scales much smaller 
than the dominant velocity length scale, as in the simulations 
of 11911 . In simulation D, P m = 1 and the field is amplified 
from an initially weak seed magnetic field by turbulent veloc- 
ities sustained by maximally helical random forcing. When 
the dynamo growth of the magnetic field saturates, the domi- 
nant magnetic field length scale is comparable to the dominant 
velocity length scale, as in the simulations of 12(11 . "Simu- 
lation" A2rp is obtained by assigning each Fourier mode in 
simulation A2 a random phase without changing the mode's 
amplitude. In simulations Al, A2, A2,p, and D we set k/Ib 
equal to the maximum of kEt,(k), where E/,(k) is the power 
spectrum of the magnetic field [the total magnetic energy is 
/ E[,(k)dk]. In simulation B we take k/Ib to be the maximum 
of kE v (k), where E v (k) is the power spectrum of the velocity 
field. We set 7t/Z ( / equal to the maximum of k i Et,(k). 
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TABLE I: 6fi/| (B) | is the ratio of the rms fluctuating field to the mean 
field, H m is the magnetic helicity divided by the maximum possible 
magnetic helicity at that level of magnetic energy, lg /lj is the ratio 
of outer scale to inner scale, and a is the coefficient in equation Q. 

We take a snapshot of the magnetic field in each simulation 
and introduce 2,000 pairs of field-line tracers whose initial 
separation ro is perpendicular to the local field. We use linear 
interpolation to obtain the magnetic field between grid points 
and employ second-order Runge-Kutta to integrate field-lines. 
To improve the statistics in simulation Al, we use 2,000 field 
lines in each of seventeen snapshots of the field separated in 
time by an interval QAIb/u, where u is the rms velocity. In 
simulations A2 and A2 lp , we use 20,000 field lines in each of 
five snapshots of the magnetic field separated in time by an 
interval 0.21b /u. We iteratively reduce the length step in the 
field-line integrations to achieve convergence. 




r/l„ 

FIG. 1 : The probability distribution function of field-line separations 
in simulation Al with ro = 2~ 9 /#. 



The probability distribution function (PDF) of field-line 
separations in simulation Al for ro = 2~ 9 Ib is plotted in fig- 
ure ^for different values of the distance z a field-line pair is 
followed along the magnetic field. [The probability is per unit 
log 10 (r/Zg), so that the probability that r lies in some interval 
is proportional to the area under the plotted curves.] The PDF 
is highly non-Gaussian, and the tail of the distribution dom- 
inates the growth of (r) when l ( [ < (r) < Is- The maximum 
kurtosis (r 4 )/(r 2 ) 2 is - 500 and occurs for z — 0.2. We find 
that decreasing ro increases the maximum kurtosis. 

The growth of (r) with z in simulation Al is shown in fig- 
ure|2]for several values of ro. There are three stages of growth 
similar to those described in p revious theoretical treatments of 
field-line divergence 1 1 CM. 1 1 111 : (1) an initial stage of exponen- 
tial growth when ro <C Id, (2) approximate power-law growth 
with r oc ■£ for (r) < l B , and (3) (r) oc z l / 2 7 for (r) > l B . How- 
ever, some aspects of our results differ from previous stud- 
ies. Within a single simulation, a increases with decreasing ro, 
which is probably related to the increasing prominence of the 
tail of the PDF as ro is decreased. In addition, stage 2 with 
(r) oc z" begins for (r) < Id, perhaps because the field-line 
pairs with largest r, which dominate the growth of (r), sat- 
isfy r > Id before (r) > lj. 

We seek to test the qualitative prediction of 
that (z s ) and Ls asymptote to a value of order a few lg as ro 
is decreased towards l d in the large-/g/^ limit. In figure [5] 
we plot (z s ) for simulation Al and simulation A2. The lower- 
resolution data of simulation Al suggest the scaling (z s ) x 
ln(Zfi/ro) for Id < ro < 0.25Zg, in contradiction to the theo- 
retical treatments. On the other hand, for simulation A2, the 
curve through the data is concave downward for lj < ro < Ib- 
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FIG. 2: Mean field-line separation (Ar) as a function of distance / 
travelled along the magnetic field in simulation A for different values 
of r . 



Moreover, figure [5] shows that the simulation Al data, and 
probably also the A2 data, have not converged to the high- 
Reynolds-number values of (z s ) for Zb/16 < ro < Ib, values 
of ro that are within the inertial ranges (Id to l B ) of both sim- 
ulations. Also, the slope d{z s ) / 't/(ln(/g / Vo)) for both ro < Id 
and ro < /b/10 decreases significantly when fa/Id is doubled. 
A comparison of the data for Al and A2 thus suggests that in 
the large-Zg/Z^ limit (z s ) asymptotes to a value of order sev- 
eral l B as ro is decreased towards Id, as in 

HE El- The 

same comments apply to the data for Ls, which are plotted in 
figure 0] 

We note that for ro = l B , (z s ) and Ls are by definition 0. 
The numerical-simulation data points that appear to be plotted 
above /g/ro = 1 actually correspond to ro just slightly smaller 
than Ib, indicating that (z s ) and Ls are discontinuous at ro = Ib 
in the numerical simulations. The reason is that for ro just 
slightly less than Ib, some fraction of the field line pairs are 
initially converging and must be followed a significant dis- 
tance before they start to diverge. 

The thermal conductivity Kt in galaxy-cluster plasmas 
scales approximately like the diffusion coefficient of ther- 
mal electrons (3 [3]. For clusters, Is/Id — h/Pi — 10 13 , 
where p,- is the proton gyroradius [21]. In the \arge.-l B /ld 
limit, the numerical simulations and theoretical models indi- 
cate that Ls asymptotes to a value of order several Is as ro is 
decreased towards Id, and Ls is not expected to increase ap- 
preciably as ro is further decreased from Id = p, to p e . Thus, 
Ls(p e ) ~ L$(ld). We take simulations Al and A2 to be our 
best models of a galaxy-cluster magnetic field. We find that 
Ls(Z d ) ~ in Al with l B /l d = 23, andL s (/ rf ) ~ 10- 
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FIG. 3: The average distance (z s ) that a field-line pair must be fol- 
lowed before separating by a distance lg as a function of initial field- 
line separation ro for simulations Al (small circles) and A2 (large 
circles). 
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FIG. 4: The dependence of Ls on tq. 



in A2 with h/ld — 50. If Ib is redefined so that 2k/ Ib (instead 
of k/Ib) equals the maximum of kEi,(k), then L$(ld) — 71b 
in Al and Ls(ld) — 6.51b in A2| 14]. Since the definition of l B 
contains an arbitrary constant of order unity, it is not clear 
which of the two definitions leads to the more accurate pre- 
diction of Kj. Given this uncertainty, extrapolating these re- 
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suits to the large-Zg /l c ( limit suggests that Ls('rf) ~ 5 — 10/g in 
intracluster plasmas. Diffusion alongthe magnetic field can 
be suppressed by magnetic mirrors 0, [Till and wave pitch- 
angle scattering 1 22], but in clusters the Coulomb mean free 
path is sufficiently short that neither of these effects is very 
important. 11 lL 1 1411 Thus, the parallel diffusion coefficient of 
thermal electrons is comparable to the thermal-electron diffu- 
sion coefficient in a non-magnetized plasma, Do. In the static- 
magnetic-field approximation, equation (|6jl thus implies that 
D/Dq ~ 0.1 — 0.2 for thermal electrons in intracluster plas- 
mas, and that Kj is reduced by a factor of ^5 — 10 relative 
to the Spitzer thermal conductivity Ks of a non-magnetized 
plasma. Heating of intracluster plasma from thermal conduc- 
tion with Kt ~ 0.1 — 0.2ks would be sufficient to balance ra- 
diative cooling in a some but not all clusters I23I E4I1 . More 
work is needed to determine the validity of the static-field ap- 
proximation in clusters, to determine the factors of order unity 
that have been neglected in our phenomenological treatment 
of thermal conduction, and to clarify the effects of turbulent 
diffusion and turbulent resistivity on heat conduction in strong 
MHD turbulence. 
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